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Cosmic Microwave Background polarization and the 
ionization history of the Universe 
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ABSTRACT 

We point out that polarization measurements as planned for the upcoming PLANCK 
mission can significantly enhance the accuracy of cosmic parameter estimation com- 
pared to the temperature anisotropy spectrum alone. In order to illustrate this, we 
consider a standard cosmological model and several modifications that adjust one pa- 
rameter each to fit the recently published Maxima- 1 data. While all models produce 
acceptable fits as far as the power spectrum is concerned, their corresponding polar- 
ization spectra differ widely. The strongest differences are expected for a model with 
delayed recombination, reflecting the fact that polarization measurements are most 
sensitive to the processes governing the epoch of recombination. 
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■ Cosmology:observations 



Observations of the Cosmic Microwave Background 
(CMB) are fundamental for our understanding of one of the 
most important epochs in the early history of the Universe. 
After a number of success ful experiments, among them 
CORE (pmnett et al. 199e|), BOOMERANG (|dc Bcrnadi s 



et al. 2001 



et al. 2|000|), Maxima- 1 ([Hanany et al. 200(]|), CBI dPadTi 



the value expected from standard Big Bang nucleosynthe- 
sis (SBBN) supported by He^ and deuterium mass frac- 
tion measurements. These discrepancies stimulate discus- 
sion about possible alternatives that either alte r SBBN, 
e.g. by non-vanishing neutrino che mical potentials (Esposito 
et al. 2001; Mangano et al. 200C), or add additional non- 



to mention just a few, much attention is now 
focused on the estimation of cosmological parameters such 
as the densities of baryons Q,b and dark matter Qm and cos- 



mologi( ]aI constant Q,/^. In addition, preliminar information 
about the Hubble constant Ho = IQQh, the power spectrum 
of initial adiabatic perturbation, and the ionisation history 
of the Universe are available. 

Fitting the CMB anisotropy power spectrum to the 
observational BOOMERANG and Maxima- 1 data confirms 
the prediction of the standard inflation scenario that the 
Universe is flat, i.e. that the total amount of matter and 
vacuum energy has the critical value - to an accuracy of 
10%, Q. b + + = 1. However, as was already pointed 



standard parameters to the cosmological model, such as a 
tilt of the power spectrum ( Tegmark fc Zaldarriaga 2000a ) 
or seco ndary ionization of hydr ogen ( [Tegmark fc Zaldarriaga 



2000b; schmalzing et al. 2000). AU these factors should be 



out by e.g^j_^c^mar^^^^aldarria^^ (2000b) ; White et 



al. (2000)1 ; [Lesgourgues fc Peloso (2000)i the data further- 
more yield a baryon fraction Q,iyh significantly larger than 



taken into account in the reconstruction of the history of 
the cosmological expansion, and therefore in the develop- 
ment of methods for cosmological parameter estimation for 
future high precision experiments like MAP and PLANCK. 

So far, however, little attention has been paid to a cru- 
cial part of the modelling of the CMB anisotropy spectrum 
- the history of hydrogen and helium recombination at red- 
shift z ~ 10'^. The standard model of hydrogen reco mbi- 
nation in the b aryonic Universe was suggested by (Pee- 



bles 1968) and ( Zel'dovich et al. 1968) and generaliz ed to 
t he dark matter mod els by (Zabotin fc Naselsky 1982) and 



(lones fc Wyse 1985). A more precise model of recombina- 
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tion taking into account ffe* was developed by (Seager et al. servational data of the Maxima- 1 experiment ( at least for 



199S). Recently, in order to test the sensitivity of methods 



and ex;)eriments to the detailed ionization history around 



^, 50) . Usi ng the offset lognor mal approach of Knox et al 

I 2 



(1998) 



recombination, this standard model of recombinat i on has 
been modif ied in various ways ( Avelino et al. 200C ; [Battyc 
et al. 2(t0lt [Hannested fc Scherrer 2001| ; [Landau et al. 20001 ). 
(Peebles et al. 200C) discuss some modification of the 



and 



Bond et al. (2000) we find a x of 20.2 for 10 



primordial ionization history of the hydrogen plasma at 
z ~ 10^ in the framework of a delayed recombination model, 
that cou ld arise in certain scena,rios with a decaying heavy 
part icle ( Sarkar fc Cooper 198£ jcott et al. 1991 Ellis et al 



degrees of freedom allowing this model to be rejected with 
more than 95% confidence. The effect of reionization is to 
suppress the primary anisotropies (the angular power spec- 
trum) at £ ^ 50 compared to COBE scales {£ ~ 10) due to 
electron scattering. A model with z^o ~ 15 and all other pa- 
rameters remaining unchanged is also shown in Figure ^ It 
is clear that this model provides a conside rably better fit to 
the dat a (x^ = 9.8) - this fact was used by Schmalzing et al] 



1992; D oroshkevich fc Naselsky 200l|) , primordial black holes (2000) to argue that under certain simplifying assumptions 



evaporation ( Naselsk ;;j^Ji)7^ Naselgk y fc Polnarev 1987 ), or 
cosmic string wakes ( |^^Ueret''aLT99£ ) . Ii i addition to the 



10 cosr nological parameters discussed by ( Tegmark fc Zal 



darriagi 2000a) in connection with the CMB anisotropy 
power spectrum, the possible delay of recombination intro- 



duces another one, namely Sa, related to the production of 
additional Lyman-a resonance photons in the background 
of primordial radiation. Their fraction Ua is determined by 
the simple equation = £aH{t)nH, where H{t) is the 
Hubble parameter and nn is the neutral hydrogen concen- 
tration. The infiuence of such additional Lyman-a photons 
on the rate of ionization is fairly straightforward. In com- 
parison to standard recombination, they compensate for the 
decrease in the number of resonance quanta due to cosmo- 
logical expansion and lead to an increase of the fraction 
of ionization, which becomes larger with the £a parameter. 
This process increases the dissipation length for initial adia- 
batic perturbations and therefore suppresses the secondary 
Doppler-peaks. It also shifts the anisotropy power spectrum 
toward lower i by increasing the acoustic horizon at last 
scattering. Together, the above-mentioned effects lead to 
significant alterations of the CMB temperature anisotrpy 
spectrum. However, an even more pronounced effect can in 
principle be seen in the CMB polarization power spectrum, 
which is more sensitive to the width of the last scattering 
surface. Therefore, measurements of the CMB polarization 
could be a critical test for investigation of the ionization 
history of the Universe. 

There is one additional reason to discuss the possible 
transformations of the CMB polarization power spectrum 
due to a more complicated ionization history of the Universe 
at the period z ~ 10^ and at the later epochs: In (Schmalz- 



ing et i\. 200C) it was discussed how a "conventional" cos- 
mological model (in particular, with Q.bh^ = 0.019 as de- 
duced from SBBN, fim — 0.3 resulting from this and the 
observed baryonic fraction in clusters of galaxies and stan- 
dard recombination) can fit the Maxima- 1 data provided 
that the universe reionized at redshifts z-^a ~ 15 — 20 corre- 
sponding to electron scattering optical depths r ~ 0.1 — 0.2. 
This effect is illustrated in Figure |l] which shows the angular 
power spectrum Ct of various models as we ll as the Maxima- 
1 data taken from Hanany et al. (2000) (the Boomerang 



and CBI data are also sho wn, however we shall in t his pa- 
per adopt the approach of Schmalzing et al. (2000} and fit 



to the Maxima- 1 data only - as was discussed in Schmalz- 
ing et c [l. (2000) very similar results are obtaine d by fitting 



to thcrtombmcd Jjoomcrang+Maxima-i data ot jattc ct al 



(2000) [)ni lij jjocn thai a CUBE-noi'mafizod, convoiilioiiaF 
model Ivith Hb/i^ = U.019, = 0.3, Oa = U.V, r = and 



other parameters given in table |l| tends to lie above the ob- 



the Maxima- 1 data can be used to place a lower limit on the 
redshift of reionization Zre > 15 (8) at the 68% (95%) confi- 
dence level. On the other hand, a z-^c much larger than 10-15 
is probably not likely, at least in the framework of current 
CDM structure formation scenarios 



see, e.g.^ 



Haiman fc 



Loeb (1998) (although reionization initially must be patchy. 



their work also indicates that the approximation of a sharp 
transition of the ionization state of the universe from al- 
most neutral to fully ionized used in our models should be 
a reasonable first approximation). 

Even keeping Q.bhP' and f2m fixed and assuming stan- 
dard recombination there are, however, (at least) two alter- 
native ways of suppressing the anisotropy power spectrum 
Ct Sit I ^ 50: either by including tensor modes (gravitational 
waves) or by assuming a "red-tilted" initial power spectrum 
(and in both cases t = Q). These models are also shown in 
Figure |l| - the additional parameter of the models has been 
determined by requiring that the amplitude of the angular 
power spectrum at the first acoustic peak (at I ^ 200) is the 
same as that of the model with reionization at z^a ~ 15, de- 
scribed above. The additional parameter is given in table ^ 
for the two models and, as can be seen from Figure |l[ both 
models result in respectable fits to the data (with x ~ 9.5 
and 6.8, respectively). 

One potential way of distinguishing between the above 
three types of well fitting models is by using observations of 
polarization anisotropies. The polarization power spectrum 
for all models discussed above is shown in Figure ^. The most 
notable difference between the reionization model and the 
others is the small "bump" at ^ 5, which is clearly seen, 
when the power-spectra are plotted on a logarithmic scale. 
However, in the displayed linear plot the bump is not seen 
and it is so small (< IfiK^) that it is unlikely that any of 
the currently planned experiments (including the PLANCK 
mission) will be able to test the existence of such a feature. 

We now turn to the delayed recombination model^ 
We adjust the free parameter Sa by normalizing to the 
Zre = 15 model at the first Doppler peak, as above. As 
seen from Figure m the model provides a good fit to the 
Maxima-1 data (x = 6.8 ) and agrees with th e recent CBI 
data point published by Padin et al. (2001) at the range 
£ ~ 1000 — 1500. In Figures |^ and |^ we show the predicted 
polarization power spectrum and temperature-polarization 
cross correlation power spectrum for the various models on 



In o rder to calculate numerical spectra for the 



Peebles et al 



(2000' mo dels, we generalize the RECFAS T code from Seagcr et 
and the CMBFAST code from 



301299) 



(1996] to incorporate the additional parameters 



Seljak & Zaldarriaga 
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model 


parameter 


value 


•^Maxima 


Standard 






20.17 


Late reionization 




15 


9.77 


Tensor modes 


Qt/Qs 


0.15 


9.52 


Tilted initial spectrum 


n 


0.95 


6.79 


Delayed recombination 




7 


6.75 



Table 1. values for comparing the Maxima data to several 
variations of a standard cosmological model. All models share 
the parameters n^/i-^ = 0.019, Q/^ = 0.7, and h = 0.65. There is 
neither a contribution from space curvature [Qfc = 0), nor from 
hot dark matter (f2j/ = 0) in any of the models. 



4000 



C^aooo 




Figure 1. The temperature power spectrum predicted for the 
various models summarized in Table |l| The lines show predictions 
obtained with a modified version of CMBFAST. The solid line 
corresponds to the standard model, while delayed recombination 
and tilted spectrum are shown in long and short dashed, respec- 
tively. Both the model with late reionization and the one includ- 
ing tensor modes shown in dotted lines, because the curves nearly 
overlap anyway. The various symbols represent data points of the 
Maxima-1 experiment (solid circles), Boomerang (open squares), 
and CBI (nothing). 



a linear scale. At large £ (~ 1300 — 1400) the difference be- 
tween the late recombination model and all the other models 
discussed becomes considerable. In fact, it is so large that 
the forthcoming PLANCK mission is expected to be able 
to discriminate between the two classes of models. Hence, 
in this case polarization measurements are expected to add 
valuable information to the one obtained from temperature 
anisotropy measurements. 




Figure 2. The po 

same as in Figure |l 



arization power spectrum. Line styles are the 



100 - 




Figure 3. The cross correlation between temperature and po- 
larization. Also here, each model is represented in the same line 
style as in Figure 0. 
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